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The highly excited stretching vibrational energy levels and the intensities of infrared 
transitions in tetrahedral molecules are studied in a U(2) algebraic model. Its applica- 
tions to silane and silicon tetrafluoride are presented with smaller standard deviations 
than those of other models. 
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I. INTRODUCTION 



In recent years, algebraic models, such as Lie algebraic methods (Iachello and Levine, 1995; Bijker 
et al., 1995) and the boson-realization model (Ma et al., 1996), have been proposed for descriptions 
of vibrations, rotations, and rotation-vibration interactions in polyatomic molecules. Lie algebraic 
methods for diatomic molecules have been modified by the corresponding quantum algebra (Alvarez 
et al., 1994; Chang, 1995), and the boson-realization model have been developed for the higher 
vibrational states of polyatomic molecules in terms of q-deformed oscillators (Xie et al., 1996; Hou et 
al, 1997). 

In Lie algebraic approaches, U(4) and U(2) algebraic model have been extensively used. The U(4) 
model took the rotation and the vibration into account simultaneously but became quite complicated 
when the number of atoms in a molecule increased to larger than four. The U(2) model was par- 
ticularly successful in explaining stretching vibrations of polyatomic molecules such as benzene-like 
and octahedral systems (Iachello and Oss, 1991; Chen et al., 1996). This model was extended to deal 
with both stretching and bending vibrations in triatomic molecules (Frank et al., 1996). Recently, a 
U(5) algebraic model was introduced for higher excited stretching modes and infrared intensities of 
tctrahcdral molecules (Leroy et al., 1996; Leroy and Boujut, 1997). However, the U(5) model was less 
feasible than the U(2) model when the bending vibrations were considered. 

In this paper, we will use the U(2) algebraic model to study the stretching vibrations and intensities 
of infrared transition of silanc S1H4 and silicon tetrafluoridc SiF4. Our results are quite good in 
comparison with those of the U(5) algebraic model for SiELt (Leroy et al., 1996) and the local mode 
model (Delia Valle, 1988). Results for stretching and bending vibrations in the U(2) model and the 
boson-realization model will be presented in a subsequent publication. 

II. U(2) ALGEBRAIC MODEL 

For a tetrahedral molecule XY 4 , we introduce four U(2) algebras to describe the vibrations of four 
X-Y bonds. The molecular dynamical group is 

[/i(2)®[/ 2 (2)®[/ 3 (2)®[/ 4 (2), 

where each Uj(2) (1 < i < 4) is generated by the operators {N i} J+^J-s, Jo,i}, satisfying the following 
commutation relations (Frank et al., 1996): 

[N i ,\ j ]=Q, [Ni, J ±J ] = 0. 

where Ni is related with the Casimir operator of U(2): 

2 + J+,iJ-,i + J-,iJ+,i = Ni{Ni/2 + 1). 
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Denote by the number of quanta in the ith bond. The local basis states for each bond are labeled 
by the eigenvalue N of Ni and v i7 and written as \Ni,Vi). Their products provide the local bases: 

\N 1 ,v 1 )\N 2 ,v 2 )\N 3 ,v 3 )\N 4 ,v 4 ) ee \N u Vi). 

where those N are equal to each other, N — N, due to equivalent bonds . 
There are three kinds of 0(2) invariant combinations of those generators: 

H l = (J+,iJ-,i + J-,i,J+,i)/2 - N/2, 

Ha = 2J ,iJ j - NiNj/2, ijtj, (1) 
Vij = J+,iJ-,j + -I i^j- 

Their matrix elements in the local bases are given by Frank et al. (1996). The operator Hi corresponds 
to the energy of the ith Morse oscillator. The operators Hij describe the anharmonic terms with the 
type ViVj, while the operators Vj describe the interbond couplings which, in configuration space, are 
of the type r^j, where and r 3 are the displacement vectors of bonds i and j from their equilibrium 
positions. 

The Hamiltonian, if restricted at the quadratic terms, is expressed in terms of those three kinds of 
operators as follows: 

H = Ai Hi + A 2 £ + A 3 % + A 4 (^) 2 
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where all A's are the coupling parameters. The Hamiltonian preserves the quantum number V = J2 v i- 
We now apply this model to study the stretching vibrational spectra of SitLi and SiF4. The calcu- 
lation for energy levels has been greatly simplified since the symmetrized bases are used. The boson 
number N is taken to be 60 for S1H4, and 100 for SiF4. 



Table I 
Table II 



For SiH4, we choose five parameters, the same as the number of parameters used by Lcroy et al. 
(1996), to calculate the vibrational levels. It was found that these five parameters Xj, 1 < j < 5, 
can give better results. That means to set Aj=0, 6 < i < 18. Fitting the observed data for S1H4 
from the compilation of Leroy et al. (1996), we obtain the parameter values and the calculated energy 
levels, listed in Table I. For SiF4, the better results can be obtained in terms of only three parameters 
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Ai, A2, and A3. Its observed data from McDowell et al. (1982), calculated values, and corresponding 
parameters arc also listed in Table I. For comparison, the levels for the two molecules calculated by 
the local model model (LMM) (Delia Valle, 1988) and the recently calculated results for SitLi by the 
U(5) model (U(5)M) (Leroy et al., 1996) are given in Table I, together with their standard deviations 
(SD). The calculated vibrational energy levels given in Table I are the differences between the observed 
values and the calculated ones. 

III. INTENSITIES OF INFRARED TRANSITION 

In the following we will introduce infrared transition operator. The calculated intensities can be 
used to check assignments and to study the intramolecular energy relaxation in tetrahcdral molecules. 

For the considered systems, the infrared active mode is F 2 . The absolute absorption intensities 
from a state v' to v are given by 

Pw = \(v\f x \v')\ 2 + \(v\f y \v')\ 2 + \(v\f z \v')\ 2 , 

where v vv i is the frequency of the observed transition, T x , T y , and T z correspond to the three com- 
ponents of the infrared transition operator T, and the state \v) denotes \Ni,Vi) for short. All other 
constants are absorbed in the normalization of the operator T. The three components of T are 

T x = a (ii - i 2 + i 3 - £ 4 ), 

T y = a {h-h-ia+h), (4) 
T z = a {h+h-i 3 -U), 

where a is the parameter. The matrix elements of ti are taken to be (Iachello and Oss, 1993) 

{N l ,v l \U\N ll v[) = exp(-0i\vi-vt\). (5) 

Those Pi for equivalent bonds are equal to each other, and denoted by a common symbol (3. 

The two parameters in the transition operator of (4) and (5) will be determined by fitting observed 
data. Due to the Wigner-Eckart theorem, it is sufficient to calculate only the z component of the 
transition operator, T z , in the symmetrized bases. In Table II we only list part of calculated intensities 
to compare with known observed data, from that the two parameters are determined. The standard 
deviation in our fitting is 1.265 for SiH 4 , and 2.512 for SiF 4 , while that was 1.415 for SiH4 by Leroy 
et al. (1996). 

IV. CONCLUSION 

We have used a U(2) algebraic model for studying the stretching vibrations and infrared intensities 
of a tetrahedral molecule. The model Hamiltonian and the model transition operator have provided 
better fits to the published experimental data of silanc and silicon tetrafluoride than the local mode 
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model (Delia Valle, 1988) and the U(5) algebraic approach (Leroy et al., 1996) to silane. In addition, 
U(2) algebraic model can be applied to the stretching and bending vibrations of other medium-size 
and large molecules. Furthermore, this model can be studied by the corresponding quantum algebra 
(Bonatsos and Daskaloyannis, 1993). Investigations on those subjects are under way. 
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Table I. Observed and calculated energy levels, and parameters (in cm x ) 





SiH 4 








SiF 4 




E fc s 


LMM 


U(5)M 


U(2)M 




LMM 


U(2)M 


2186.8730 


0.5339 


0.7378 


0.4476 


800.6 


-5.1243 


-1.5226 


2189.1901 


-0.0652 


-0.5527 


-0.0765 


1031.3968 


-3.1250 


1.7942 


4374.5600 


0.3247 


-0.1007 


-0.0503 


2059.1 


-0.2942 


4.4876 


4308.3800 


-1.4258 


-0.8963 


-0.6650 


1828.17 


-1.5378 


-0.0686 


4380.2800 


1.7690 


0.2231 


-0.4878 


2623.8 


4.2976 


0.9521 


4378.4000 


1.3154 


0.1421 


-0.6465 


3068.5 


0.5720 


-4.1160 


4309.3485 


-0.5244 


0.0617 


0.2095 


SD 


3.090 


2.066 


6496.1300 


0.7661 


-0.0028 


0.0190 








6361.9800 


-2.1182 


0.1450 


0.5381 








6500.3000 


2.5976 


1.4384 


1.3330 








6502.8800 


2.2612 


0.4108 


0.4409 


Par am. 


SiH 4 


SiF 4 


6362.0800 


-2.0186 


0.2448 


0.6372 


Ai 


30.2883 


1.1974 


6497.4810 


0.9458 


-0.4104 


-0.0605 


A 2 


-2.2592 


-2.8473 


6500.6000 


1.0278 


-0.1297 


-0.05781 


A 3 


-0.0147 


-0.5687 


8347.4000 


-3.8164 


0.1031 


0.2949 


A 4 


-0.0021 




10267.2200 


-3.9023 


0.4100 


0.1748 


A 5 


-0.0003 




12121.2000 


-2.7061 


-0.6399 


-0.1973 








13914.4000 


4.8506 


0.2227 


0.3037 








SD 


2.465 


0.610 


0.573 









TABLE II. Observed and Calculated Relative Intensities 





a 


8.130 


SiH 4 


(3 


0.866 


SiF 4 

a 37.977 (3 3.540 


Obs. a 
Calc. 6 
Calc. 


10267.04 12121.40 13914.10 15645.57 17318.60 
100 21 2.4 0.6 
530 100 20.18 4.2 0.9 
479 100 20.33 4.04 0.79 


1029.60 1828.24 2054.60 2622.84 3072.62 
5000 7 1.2 0.015 0.015 

5000 3.79 4.14 0.005 0.003 



a Observed data for SiH 4 from Bernhcim et al. (1984), and for SiF 4 from McDowell et al. (1982). ^Calculated by Leroy 
et al. (1996). 
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